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We introduce a class ol exactly solvable SO(n) symmetric Hamiltonians with matrix product 
ground states. For an odd n > 3 case, the ground state is a translational invariant Haldane gap 
spin liquid state; while for an even n > 4 case, the ground state is a spontaneously dimerized state 
with twofold degeneracy. In the matrix product ground states for both cases, we identify a hidden 
antiferromagnetic order, which is characterized by nonlocal string order parameters. The ground- 
state phase diagram of a generalized SO(n) symmetric bilinear-biquadratic model is discussed. 
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I. INTRODUCTION 



One-dimensional quantum Heisenberg antiferromag- 
nets have a long history and show many fascinating prop- 
erties. Since the Mermin- Wagner- Coleman theorem 1,2 
forbids a continuous symmetry breaking in one dimen- 
sion, no classical Neel order can survive, even in zero 
temperature. The rigorous solutions on particular mod- 
els provide essential insights to understand the proper- 
ties of these quantum spin liquid states. For instance, 
the spin- 1/2 antiferromagnetic Heisenberg chain has a 
Bethe-ansatz solution; 3 -^ which yields a unique spin sin- 
glet ground state, gapless spin- 1/2 excitations and power- 
law decay spin correlations. Meanwhile, additional next- 
nearest-neighbor interactions can frustrate the nearest- 
neighbor antiferromagnetic correlations. The Majumdar- 
Ghosh model^ is such an exactly solvable example, which 
has a twofold-degenerate dimerized ground state, a finite- 
energy gap, and extremely short spin correlations. 

Toward the quantum integer-spin models, Haldane 
gave a striking prediction that an excitation gap occurs 
between the ground state and the excited states^ Al- 
though Haldane's argument is based on a semiclassical 
large- «S expansion, it was later verified by numerical stud- 
ies for lower-5 casesi^jS Remarkably, Affleck, Kennedy, 
Lieb, and Tasaki (AKLT) found a family of integer-spin 
chain Hamiltonians with exact massive ground states, 
which are called valence bond solid (VBS) states^ The 
VBS states preserve spin rotational symmetry, and ex- 
hibit exponentially decay spin correlations and gapped 
excitations, thus, share the key features of Haldane gap 
spin liquid states for the quantum integer-spin Heisen- 
berg antifcrromagnets. Although no true long range or- 
der exists, den Nijs and Rommelseii observed a hid- 
den antiferromagnetic order in the 5=1 VBS state, 
and introduced a set of nonlocal string order parame- 
ters to provide a faithful quantification of the S = 1 
Haldane phase. The string order in S = 1 VBS state, 



Haldane gap, and the fourfold degeneracy in an open 
chain can be understood by a hidden Z2 x Z2 symme- 
try breakingi 12 ' 13 i 14 However, a nonlocal string order pa- 
rameter that reflects correctly the hidden Z$+i x Z$+i 
symmetry of the higher- S VBS states remains an open 
problem^ii^iie 

Beside the studies on S'C/(2)-symmetric spin chains, 
quantum spin systems with higher symmetry also at- 
tract much attention. For instance, the Bethe-ansatz 
method for SU(2) Heisenberg chains can be generalized 
to models with SU (n) symmetry^ It has been argued 
that such an SU (4)-symmetric model can be achieved 
in electronic systems with two-fold orbital degeneracy 
at quarter filling.— Meanwhile, Affleck et al*& first dis- 
cussed the extension of VBS states to £J7(2n)-mvariant 
extended VBS states, which break lattice translational 
symmetry and charge-conjugation symmetry but remain 
invariant under the combined operation of these two sym- 
metries. Furthermore, Greiter et al£Q studied the SU(n) 
spin chains with exact valence bond solid ground states. 
Along with the rise of cold atomic physics in optical lat- 
tices, Chen et alS^ constructed an SU(4) Majumdar- 
Ghosh model with exact plaquette ground states by us- 
ing spin-3/2 fermions. Very recently, Arovas=* explored 
a family of novel SU(n) simplex solid states, which are 
natural generalizations of SU(2) VBS states of AKLT 
models. Besides these models with SU(n) symmetry, 
Schuricht and Rachel 23 considered Sp(2n) VBS states 
and their parent Hamiltonians. 

In this paper, we will introduce a class of SO(n)- 
symmetric Hamiltonians with matrix product states as 
their exact ground states. However, these SO{n) sym- 
metric spin chains show a different even-odd effect. For 
an odd n — 21 + 1 , a periodic chain has a unique ground 
state. All these SO(2l + 1) matrix product states have 
a hidden antiferromagnetic order, which is character- 
ized by string order parameters. The nonlocal unitary 
transformations are designed to explicitly reveal a hid- 
den [Z2 x Z2) 1 symmetry. The breaking of this sym- 
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metry is responsible for the Haldane gap, nonvanishing 
string order parameters, and 4 z -fold degeneracy in an 
open chain. However, for an even n = 21, a periodic 
chain has a twofold dimcrized ground state, which breaks 
translational symmetry. Nevertheless, these SO(2l) ma- 
trix product states also contain a hidden antiferromag- 
netic order. Finally, the ground-state phase diagram of a 
generalized SO(n) symmetric bilinear-biquadratic model 
is obtained. 

This paper is organized as follows. In Sec. II, the 
SO(n) algebra and the exactly solvable SO(n) symmetric 
models will be introduced. In Sec. Ill, the exact matrix 
product ground state of the SO(n) model with n — 21 + 1 
will be studied, in particular, with the examples of n — 3 
and 5. The hidden order in all these SO(2l + 1) matrix 
product states and the corresponding hidden symmetry 
are identified. Section IV is devoted to an analysis of the 
SO(n) model with n = 21, which has a dimerized ground 
state breaking lattice translational symmetry. In Sec. 
V, a generalized SO(n) symmetric bilinear-biquadratic 
model is introduced and their ground-state properties are 
discussed in detail. A conclusion is presented in Sec. VI. 



II. MODEL HAMILTONIAN 

Let us begin with a one-dimensional SO(n) chain of 
N lattice sites (N even). On each site, the local Hilbert 
space C" contains n states \n a ) (1 < a < n), which 
can be rotated within the SO(n) space via the following 
vector relations: 

L ah \n c )=i5 bc \n a )~i5ac\n h ) 1 (1) 

where L ab (a < b) are the n(n — l)/2 generators of the 
SO(n) Lie algebra. The vector relations constitute the 
n-dimensional representation of SO(n) algebra and the 
following commutation relations holdi^i 

[L a \ L cd ] = i (5 ad L bc + 5 bc L ad - S ac L bd - S bd L ac ) . (2) 

According to the Lie algebra, the tensor product of 
two SO(n) vectors can be decomposed as a direct sum of 
an SO{n) singlet with a dimension 1, an antisymmetric 
SO(n) tensor with a dimension n{n — l)/2, and a sym- 
metric SO(n) tensor with a dimension (n + 2)(n — l)/2, 
i.e., 

n®n = l® n{n- l)/2 ® (n + 2){n - l)/2 , (3) 

where the number above each underline is the dimen- 
sion of the corresponding irreducible representation. For 
n = 3, we recover the well-known Clebsch-Gordan de- 
composition 3®3 = !©3©5of two spin-1 representa- 
tions. According to the SO{n) decomposition scheme 
(|3"|), the wave functions in each irreducible represen- 
tation channel can be obtained explicitly. The maxi- 
mally entangled SO(n) singlet wave function is written 



as 77^ ELi l n<l )j' an d the wave functions of the 
antisymmetric channel are expressed as 

Finally, the symmetric channel contains n(n— 1)/2 states 
with the wave functions 
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and the rest of n — 1 states with the wave functions 

^ (V), K), - \n b \ \n») .) . (6) 

For the three SO(n) channels given in Eq. ([3]), the 
bond Casimir charge ~}2 a<b {L^ h + ^fY f° r two adja- 
cent sites takes the values 0, 2n — 4, and 2n, respec- 
tively. Together with the single-site Casimir charge 
J2 a <b(^i b ) 2 = 71 ~ 1' 0110 can wr ite the SO(n) symmetric 
bilinear interaction term as a polynomial of bond projec- 
tion operators, 

J2 L t L f = (l-n)Pi(U)- P ra(n _ 1)/2 (U) 

a<b 

+ V{n+2)(n-l)/2 {hj), (?) 

where the bond projectors Vi(i,j), V n i n -x)/2{i> j), an d 
T~'(n+2)(n-i)/2(h j) project the states of two adjacent sites 
i and j onto the three SO(n) channels in Eq. ([3|), re- 
spectively. Using the property of projection operators, 
we square Eq. (JT)) and obtain the SO(n) symmetric bi- 
quadratic interaction term as 




= (n - lfVljiJ) + Vn(n-l)/2(i,j) 



+' P(n+2)(n-l)/2 {i,j)- (8) 

Combined with the completeness relation of the projec- 
tors, 

~Pl(h j) + ' Pn{n-l)/2 (i, j) +' P(n+2)(n-l)/2 (hj) = 1, (9) 

we can express the bond projection operators with the 
SO(n) generators as 

W+^n-X)^ )) V *g \ & / 

x I Z a<b LfLf 1(10) 

V(E a <6^f)V 

Now we define our model Hamiltonian as 

HsO(n) = ^2 ^("+2)(n-l)/2 (^ i + 1), (11) 
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which is a bilinear-biquadratic Hamiltonian in terms of 
the SO(n) generators according to Eq. (TIT)]) . This model 
has exact matrix product ground states, which will be 
extensively studied below. Although the exact excited 
states are not known, we argue that there is a finite en- 
ergy gap above the ground states. For a projector Hamil- 
tonian such as Eq. (jlip , this argument can be proved rig- 
orously using a method proposed by Knabe^ who found 
that the lower bounds of energy gaps of infinite systems 
can be obtained by diagonalizing finite-size systems. 

In Sees. Ill and IV, we will discuss odd n = 21 + 1 
and even n — 21 cases separately because the nature of 
the matrix product ground state depends on the parity of 
n. Mathematically speaking, the SO{2l + 1) and SO (21) 
algebras are quite different. According to the Cartan 
classification scheme^! the SO (21 + 1) algebra belongs 
to Bi type, while the SO(2l) algebra are Di type. 



III. ODD-n CASE 

Let us assume n = 21 + 1, where I is an integer (I > 1). 
To achieve the exact ground state of model Hamiltonian 
([TTjl . one has to resort to a fascinating property of SO(n) 
algebra — the spinor representation. An elegant way to 
construct the spinor representation of SO(2l + 1) alge- 
bra is to introduce (21 + 1) gamma matrices satisfying 
the Clifford algebra {r a ,r b } = 25 a b- Then an irreducible 
spinor representation of SO (21 + 1) is immediately con- 
structed by r Qb = [r a ,r b ]/2i. The product of T a and T b 
can be expressed as r a T b = S a b + «r ab . For each lattice 
site j, if the following matrix state is introduced: 



fJ:> 



21 + 1 

E r ° 

a = l 



(12) 



then the bond product of gj at any two neighboring sites 
9j9j+l is given by 

g m+x = *E rab (K),l"Vi-|A>Vi) 

a<b 

+£k»Vi> ( 13 ) 



where the first two terms belongs to the antisymmetric 
channel in Eq. Q and the latter term is the 50(2/ + 1) 
singlet. Since Vi n+ 2){n-x)J2(j-,j + 1) performs the pro- 
jection onto the states of the symmetric (n + 2)(n — l)/2 
channel, the matrix product state defined by 

I*) = Tr (5i52 ■ ■ -9n) 

Tr(r ai r a2 



£ 

ai . . .ajv 



r ajv )|n ai n a2 ...7i ajv )(14) 



is always the zero energy ground state of the Hamiltonian 
(llip in a periodic boundary condition. This state pre- 
serves SO(2l+l) symmetry and lattice translational sym- 
metry. For an open chain, there are totally 4'-degenerate 



ground states, which can be distinguished by their edge 
states. 

To compute the correlation functions in the matrix 
product ground state, we set up a transfer matrix 
method^ by introducing 



G A = ^(n a |A|n b )(f Q (g)r b ), 



(15) 



where A is an operator acting on a single site and f a 
denotes the complex conjugate of T a . Specifically, the 
transfer matrix G is written as G = ^ Q f a <§5 T a . Then 
a two-point correlation function in thermodynamic limit 
can be written as 



(LfLf) = lim 

x 31 N- — 



Tr[(G) N ^+*~ 1 Gab 1 Gab] 



Tr(G) 



N 



(16) 



where G ab = i(T a (g> T b - T b (g> T a ). In a long distant 
limit, the two-point correlation functions of SO (21 + 1) 
generators decay exponentially as 



(LfLf) ~ exp 



\3~i\ 



(17) 



21 + 1 
21-3 



with the correlation length £ = 1/ In 

Here we note that the SO (21 + 1) symmetric model has 
a deep relation with the quantum integer-spin chains. On 
each lattice site, the (21 + 1) vectors of SO(2l + 1) can 
be constructed from the S = I quantum spin states. In 
the SU(2) spin language, the last two channels in Eq. 
© for odd n = 21 + 1 correspond to the total bond spin 
S = 1, 3, . . . , 21 — 1 and S — 2,4, . . . , 21 states, respec- 
tively. In other words, the SO(2l + 1) bond projection 
operators can be expressed using the spin projection op- 
erators Ps= m as 

P 2 p+t(»,j) - E p s=2 m -i(i,j), (18) 

m— 1 
/ 

Pag+gifoi) = E p s= 2m (*,j)- (19) 

m—l 

Thus, the role of V2i 2 +3i(hj) is to project onto nonzero 
even total spin states. Based on this property, we can 
further show that the matrix product wave function (|14p 
is also the ground state of the following quantum integer- 
spin Hamiltonian: 

H SU (2) = E E J ™ P S=2m(M + 1), (20) 
i m—l 

with all J m > 0. This model can be written as a poly- 
nomial of nearest neighbor spin exchange interactions 
Si ■ Si + i up to 21 powers and is therefore SU (2)-invariant. 
However, the ground state (|14p possesses an emergent 
SO(2l + 1) symmetry. 
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It is interesting to compare Hsu{2) with the AKLT 
model of valence bond solid proposed by Affleck et 



21 



^2 K mPs=m(i,i + l-), 
m=l+l 



(21) 



with all K m > 0. The ground state of TJaklt is also a 
matrix product state similar to Eq. (fl4|) . but the local 
g matrix for AKLT model is now a (5 + 1) X (5 + 1) = 
(7 + 1) x (7 + 1) matrixJ^ When I = 1, both models Hso(n) 
and H SU (2) become exactly the same as the 5=1 AKLT 
model if aklt i whose ground state is the celebrated 5 = 
1 VBS state. When I > 1, we emphasis that Hsu(2) 
and -Haklt differ from each other. In Sees. Ill B and 
III C, we will show that their matrix product ground 
states have very different hidden structures and belong 
to different topological phases, although they both belong 
to the Haldane liquid states. 



A. SO (3) matrix product state: S = 1 VBS 

In order to investigate the property of the SO(2l + 
1) matrix product state, we briefly review the 50(3)- 
symmetric 5 = 1 VBS state as a warm up. In this case, 
the SO (3) vectors can be represented by the 5=1 spin 
states, 

I" 1 ) = ^=(|-1>-|1», W) = -5=(|1>+|-1», |n 3 ) = |0) . 

(22) 

and the 50(3) generators are defined by spin-1 operators 
as L 12 = -S z , L 13 = S y , and L 23 = -S x . Moreover, 
the Clifford algebra is satisfied by the Pauli matrices as 
{<j a ,<j b } = 2Sab- According to Eq. (12]), the local g 
matrix can be written as 



9j 



-V2\l) 



V21-1) 



(23) 



which generates the matrix form of the 5=1 VBS state. 
Although the two-point spin-correlation functions in this 
state decay exponentially as shown in Eq. (|17[) . it has 
been observed that the upper and down spins lie alter- 
nately along the lattice, sandwiched by arbitrary number 
of non-polarized spin states. This hidden diluted anti- 
ferromagnetic order can be characterized by a nonlocal 
string order parameter first proposed by den Nijs and 
Rommelseiii 



= lim (5f TT exp(iirSP)S?) 

|j-»|-KX> 



4 

9' 



(24) 



where jj, = x, y, or z. 

On the other hand, the 5=1 VBS states on a finite 
open chain have two nearly free 5 = 1/2 edge degrees 
of freedom at the end of the chain and are thus fourfold 
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FIG. 1: Kennedy- Tasaki unitary transformation denned in 
Eq. (|25p for a typical configuration of the spin-1 VBS states. 
The hidden antiferromagnetic order is transformed to a dilute 
ferromagnetic order. 



degenerate. Both the hidden string order and the de- 
generacy in an open chain can be understood as natural 
consequences of a hidden Zi x Zi symmetry breaking. 
To manifest the hidden symmetry, the key is a nonlocal 
unitary transformation defined b y 12 ^ 13 



(25) 



j<i 



In the standard S z representation, exp (iirS x ) flips \m) 
to \—m) (m = ±1,0) and multiplies the state with a 
phase factor (—1). The physical meaning of the unitary 
transformation U can be explained as follows. For a given 
spin configuration on a finite open chain, all |0) are left 
alone and we look for the non-zero spins from the left 
to the right. Suppose there is a non-zero spin at site i, 
we count the number of |1) and | — 1) on the sites to the 
left of site i. If the number is even, we left the spin at 
site i unchanged. If the number is odd, we flip the i-site 
spin. Finally, an additional phase factor (—1) may be 
taken into account, depending on the total site number 
and each spin configuration. An example of the unitary 
transformation U on a typical configuration of the spin-1 
VBS state is shown in Fig. [TJ 

When applying the Kennedy- Tasaki unitary transfor- 
mation to the Hamiltonian, the 50(3)-symmetric AKLT 
model is transformed to a model with a discrete x 
Zi symmetry , 12 ! 13 and two of the den Nijs-Rommelse 
string order parameters become the usual two-point spin- 
correlation functions. Thus, the nonvanishing string or- 
der parameter measures the hidden Zi x Zi symmetry 
breaking of the original model. The breaking of this hid- 
den discrete symmetry leads to the opening of the Hal- 
dane gap, the hidden antiferromagnetic order, and the 
fourfold degeneracy in an open chain and thus provides a 
unified explanation of the exotic features in 5 = 1 VBS 
states. 



B. SO (5) matrix product state: A projected 
valence bond solid 

The next example is the 50(5)-symmetric matrix 
product state with I — 2. Actually, Scalapino et al. 28 
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proposed this state to describe the SO (5) "superspin" 
phase on a ladder system of interacting electrons. Here, 
it is convenient to introduce the 5*0(5) vectors by means 
of the S — 2 states, 



-J=(|-2)-|2))> 2 H-l(|2) 
^(|-1)-|D),I^) = A (|1) 



1-2)), 
1-1)), 



\n 5 } = |0) 



Moreover, we define the 50(5) gamma matrices as 



r 1 

p4 



r 2 

»a 2 , r 5 



„3 



(26) 



(27) 



Then the local g matrix can be written as 



fJj 



( \Q) j y/2\-\) i y/2\-2) s \ 

-V2|l) 3 - 10)^- V^|-2) 3 - 

V2|2) 3 -\0) j -yft\-l) 

V5|2), v^ll), |0). 



(28) 

In fact, the 50(5) matrix product state can be in- 
terpreted as a projected 50(5) VBS state (Fig. [5]). By 
using two spin-3/2 fermions, the spin-2 states can be con- 
structed as 2 ^ 

|2) = ^1 1>\ |vac) , |-2) = Vl x Vl 3 |vac) , 

2 2 2 2 

|1) = 4V>li|vac), |-l) = ^UI 4 |vac), 



|0) = -^(^t Via +^t^x)|vac) 

y2 2 2 2 2 



(29) 



where ^ creates a fermion with spin components a = 
±3/2, ±1/2. Because only site quintet (S = 2) and site 
singlet (5 = 0) are allowed for two spin-3/2 fermions on a 
single site, an extra projection has to be implemented to 
remove the site-singlet state. Owing to SO (5) ~ Sp(4), 



there exists an antisymmetric matrix 1Z 
the following properties: 



er 2 with 



n 2 = -l, v) = n- 1 =K T = -n, 



KT a K- L = {T a ) r , TZT ab n^ 1 = -(T a 



ib\T 



(30) 



Using the 1Z matrix, the 50(5) matrix product state in 
a periodic chain can be written in a projected VBS wave 
function as 

l*SO(6)> = n Vs = 2 (Nl2^L n <x^+l,0) l VaC > » 

j a/3 

(31) 

where Vs=2(j) is the site-quintet projector and 
V'] cJlapipj+i s is an 50(5)-invariant valence bond 
singlet. For an open boundary condition, the chain is 
ended with two nearly free spin-3/2 degrees of freedom 
leading to 16 degenerate ground states. Here we recall 
that the edge states of the S = 2 VBS states of the AKLT 



FIG. 2: (Color online) The schematic of a projected SO(5) 
VBS state. Each dot denotes a spin-3/2 fermion. The solid 
lines represent SO(5) singlet valence bond, and the dashed 
circles indicate the projection of two spin-3/2 fermions to form 
spin-2 site-quintet states. 



model are spin-1 degrees of freedom, which are sharply 
different from our 50(5) matrix product states. 

Similar to the spin-1 VBS state, the 50(5) matrix 
product states have an interesting hidden string order. 
Since the 50(5) algebra is rank 2, one can classify the 
states by using two quantum numbers (weights) corre- 
sponding to the mutual commuting Cartan generators 



L 12 and L 



34 



a.s 



L \mi,m 2 ) = mi|mi,TO2), 



L 34 \mi,m 2 ) = 7712 I?""!, 771 2 ) 



(32) 



These states, characterized by the SO (5) weights, are 
related to those states denoted by the usual S z quantum 
numbers as follows: 

|1,0) = |2), |-1,0) = |-2), |0,0) = |0>, 
|0,-1) = |1), |0,1) = |-1). (33) 

When we define 

r ± = l(r 3 ±ir 4 ), r ± = i (r 2 ± zr 1 ) , (34) 

the local g matrix in Eq. (|28p can be rewritten as 

9j = v^r+io,^.- ^-10,-1). + v^r+ 1-1,0^ 
+V2r_|i,o) J + r 5 |o,o) J . (35) 

By considering the property of Clifford algebra, it can be 
found that |1,0) and | — 1,0) must appear alternately in 
the 50(5) matrix product states despite arbitrary num- 
bers of |0, 0) and |0, ±1) between them. At the same time, 
|0, 1) and 1 0, — 1) also appear alternately with arbitrary 
numbers of |0,0) and |±1,0) between them. For exam- 
ple, a typical configuration of the 50(5) matrix product 
state is 

rnx : ■•• t | | | | | t •- 
777 2 : •••TOjOOOTIOOOTOjO--- 

where (T,0, |) represent \m) = (|1), |0), |-1)). This dilute 
antiferromagnetic order is in analogy with the spin-1 va- 
lence bond solid (VBS) state in terms of the S z quantum 
number, but here two quantum numbers are associated 
with the Cartan generators L 12 and L 34 . However, such 
an intriguing feature is not enjoyed by the 5 = 2 VBS 
states of Affleck, Kennedy, Lieb, and Tasaki (AKLT) 
model. Actually, the characterization scheme of the VBS 
states for 5 > 2 remains a challenging open problem. 
However, the hidden order of all 50(2/ + 1) matrix prod- 
uct states can be fully identified in a systematic and com- 
pact form. 
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C. Hidden order in the 50(2/ + 1) matrix product 
state 

Now we are in a position to identify the hidden order 
in all the SO(2l + 1) matrix product state HH), which 
is inspired from the analysis of 50(5) matrix product 
state. Since S0(2l + 1) is a rank-? algebra, one can al- 
ways choose the mutually commuting Cartan generators 
as {L 12 , L 34 , . . . , L 2l ~ 121 }. At each site, the quantum 
states are classified by the eigenvalues of these Cartan 
generators as 

L 2a - X ' 2a \m a ) = m a \m a ), (m Q = 0,±l). (36) 

Thus, the single-site states are associated with I quan- 
tum numbers {mi,-- - , mi} and they are subjected to 
the constraint 

m a mp = 0, (a^ 0). (37) 

According to Eq. (fTJ) , all the Cartan generators annihi- 
late the "extra dimension" vector |?i 2i+1 ) = |0, 0, . . . , 0). 
The other basis states can be chosen as 

|0 . . . ,m a = ±1, . . . 0) = (\n 2a ) ± i In 2 "" 1 )) . (38) 

From the property of the Clifford algebra, the hidden 
antiferromagnetic order of the ground state \^>) can now 
be identified. In any of the m a (a = 1 ~ /) channel, it 
can be shown that |m Q ) is diluted antiferromagnetically 
ordered, the same as for the 5 = 1 VBS state. Namely, 
the states of m a = 1 and — 1 will alternate in space if all 
the m a = states between them are ignored. 

This hidden antiferromagnetic order can also be char- 
acterized by nonlocal string order parameters. Similar to 
the I = 1 case, the string order parameters can be defined 
as 

j-i 

O ab = lim (L? TTcxp(iTrLf)Lf). (39) 

U-tl-oo f*. 

Since the ground state is SO (21 + 1) rotationally invari- 
ant, the above nonlocal order parameters should all be 
equal to each other. Thus, to determine the value of 
these parameters, only O 12 needs to be evaluated. One 
can compute the value of these string order parameters 
by the transfer-matrix techniques but there is an alter- 
nate intuitive approach. In the L 12 channel, the role of 
the phase factor in Eq. (|3"9")) is to correlate the finite spin- 
polarized states in the mi channel at the two ends of the 
string. If nonzero mi takes the same value at the two 
ends, then the phase factor is equal to 1. On the other 
hand, if nonzero mi takes two different values at the two 
ends, then the phase factor is equal to — 1. Thus, the 
value of O 12 is determined purely by the probability of 
mi = ±1 appearing at the two ends of the string. It is 
straightforward to show that the probability of the states 
mi = ±1 appearing at one lattice site is 2/(21 + 1) and 
thus O 12 = 4/(2Z + l) 2 . 



mi otoo|tooojto|ot 
m 2 | o 1 o o o f | o o o f o | o 

ofoottooottotof 

t0|000f|000f0|0 

ofoottooottofot 
tofoooftooofoto 

FIG. 3: (Color online) Changes of a typical configuration of 
the SO (5) ground state under the unitary transformation de- 
fined by Eq. (|41[) . Ui and U2 transform successively all mi 
and m2 states to two diluted ferromagnetic configurations, 
respectively. 

In the SO(2l + 1) Lie algebra, (L 2 " -1 - 2 ", £,2a-i,2i+i j 
and £ 2Q,2i+1 ) span an 50(3) sub-algebra in which 
exp(i7rL 2Q,2 ' +1 ) plays the role of flipping the quantum 
number m a . This exponential operator can flip the quan- 
tum numbers of m a without disturbing the quantum 
states in all other channels. This indicates that if we 
take the following nonlocal unitary transformation in the 
m a channel: 

c/ Q =n cx p(^r i,2 % 2Q,2i+1 )' ( 4 °) 

j<i 

then all the configurations in this channel will be ferro- 
magnetically ordered. Furthermore, by performing this 
nonlocal transformation successively in all the channels, 

l 

U=Y[U a , (41) 

a=l 

then all the configurations of the ground state will be- 
come ferromagnetically ordered. As an example, Fig. [3] 
shows how a typical configuration of the 50(5) matrix 
product state is successively changed under this nonlocal 
unitary transformation. 

By applying the unitary transformation (|41|) to the 
Cartan generators, it can be shown that 

i-l 

ULfU- 1 = Lf exp(z7r L f)- (42) 

i=i 

Substituting this formula to Eq. (|3"9")) . we find that 

O ab = lim (LfLf) (43) 

Thus, the nonlocal string order parameters O ab for Car- 
tan generators become the ordinary two-point correlation 
functions of local operators after the unitary transforma- 
tion. 

Under the above transformation, the SO(2l + 1) sym- 
metry of the original Hamiltonian is reduced and deter- 
mined by the symmetry of the unitary transformation 
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operators. In the m a channel, it can be shown that the 
unitary operator U a possesses only a Z 2 x Z 2 symme- 
try. Therefore, the Hamiltonian after the transformation 
has a (Z 2 x Z 2 ) 1 symmetry. This is the hidden topolog- 
ical symmetry of the Hamiltonian, associated with the 
hidden order of the original matrix product state 
When it is applied to an open chain system, the hidden 
[Z<x x Z2) 1 topological symmetry of the Hamiltonian will 
be further broken, yielding 2'-free edge states at each 
end of the chain. Therefore, the open chain has totally 
4 ; -degenerate ground states, which can be distinguished 
by their edge states. 



IV. EVEN-n CASE 

Let us assume n = 21 (/ > 2). Using the (22 + 1) gamma 
matrices, the spinor representation of the SO(2l) algebra 
can be constructed by leaving out T 2l+1 . However, we 
note that the resulting 2 i -dimensional spinor representa- 
tion generated by r° (a — 1 ~ 21) is reducible, in contrast 
to the 50(2/ + 1) algebra^ Since r 2i+1 commutes with 
all the SO{2l) generators r ah , one can construct the fol- 
lowing projection operators onto two different invariant 
subspaces: 

p± = i(i±r 2n+1 ). (44) 

For each lattice site j, we introduce the local g matrix 

as 

21 

<?.,• = £r»K>., (45) 
0=1 

then the exact matrix product ground states of the 
Hamiltonian (lll|) for n = 21 are given by 

|*±) = Tr(P ±gi g 2 ...g N ) 

= Tr(P ± r ai ...r Qjv )|n Ql ...n a ").(46) 

ai . . .ajy 

Due to the equation P±T a = T a P T (a = 1 ~ 21), we can 
observe that the states |\P±) are dimerized states and are 
connected to each other by translating one lattice site. 
Thus these two states break translational symmetry while 
they preserve the 50(2/) rotational symmetry. For an 
open chain, the matrix product ground states are 2 2 ' -1 - 
fold degenerate when combining dimerization and edge 
states. 

The static correlation functions can be computed by 
the transfer-matrix method as well. We find that the 
50(4) matrix product states have only nearest-neighbor 
correlations (L^Lf^) = —1/4 and the correlation 
length is zero. For / > 3, the two-point correlation 
function (Lf b Lj h ) has an exponential tail at a large 
distance, as in Eq. I|17p. and the correlation length is 



Although these two-point correlation functions of the 
S0(2n) matrix product states are short range, there 
is a hidden antiferromagnetic order, similar to the 
50(2/ + 1) matrix product states. Because SO(2l) is 
a rank-/ algebra, the Cartan generators can be chosen 
as {L 12 , L u , £20-1,2/ j. Thus ^ the stateg can be cliar _ 

acterized by the 50(2/) weight using Eq. (138|) . In the 
50(2/) case, the only difference is the absence of the ex- 
tra dimension vector |0, 0, . . . , 0) annihilated by all Car- 
tan generators. To measure this hidden order, one can 
use the string order parameter in Eq. ([55)) . A straight- 
forward calculation shows the value of these string order 
parameters that are given by 

1 

O ab = lim (Lf l[exp(mLf)Lf) = n . (47) 

|j-j|^oo - LJ - J L z 

r—i 

Here we note that translational symmetry breaking 
distinguishes the 50(2/) case from the 50(2/ + 1) case, 
while the latter belongs to the Haldane spin liquid class. 
This is an interesting even-odd effect. Furthermore, one 
may expect that their low-lying excitations are also very 
different. The low-energy excitations in the 50(2/ + 1) 
Haldane liquid are magnons, while the 50(2/) systems 
are soliton-like excitations connecting the two dimerized 
states. Although the exact results of the low-lying excita- 
tions do not exist, more evidence comes from the SO (4) 
case. 



A. SO (4) matrix product state: A staggered 
spin-orbital crystal 

The 50 (4) case is somewhat special because it can 
be factorized as 50(4) ~ 5/7(2) x 5/7(2) with C 2 <g> C 2 
being its vector representation. Namely, one can consider 
a spin-orbital 5 = T = 1/2 coupled chain or equivalently 
a spin-1/2 two-leg spin ladder to implement the SO (4) 
vectors and generators in Eq. (fTJ). We find that it is 
convenient to introduce the four vector states as 

±itt/4 

K 2 ) = ^-(iu) + iu», 

K' 4 ) = ^^(||,T) + |T,|», (48) 

where the first index in \<t,t) denotes the spin direction 
while the second one is the orbital direction. Moreover, 
the 50 (4) generators are defined by 

L 12 = _ T z _ L 13 =T x _ gx^ L 14 = _ T y _ gy^ 

L 23 = T v_ S v j L 2i = T x + S x , L 34 = T Z -5 Z ,(49) 

where 5 M and T M (/1 = x, y, or z) denote the spin and 
orbital degrees of freedom, respectively. Alternately, 5 M 
and T M can be viewed as the spin operators in the upper 
and lower leg of a two-chain ladder. 
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(a) 



(b) 



FIG. 4: (Color online) The schematic of the twofold- 
degenerate staggered spin-orbital crystal states (a) |\&+) and 
(b) |*_). 



A convenient choice of T matrices for 5*0(4) spinor 
representation is given by 



r 1 = ^ 

r 3 = a - 



cr 3 , T 2 = -a 1 (g) a 3 , 
a 2 , T* = -a°»a 1 . 



(50) 



The invariant subspace projector is P± = (l±T 5 )/2 and 
r 5 = a 3 ® a 3 . A little calculation shows that the local g 
matrix is given by 







o 



9 j = 



-IT, I), o 
i|T,T> 3 - o 
o — < IT, T>^ 



o 



iu>, 



-i\i,i) 3 
-\U)j 
/ 



(51) 



up to an unimportant normalization factor. 

In this case, the twofold-degenerate ground states have 
an intuitive meaning, which becomes clear when the local 
Hilbert space is represented by a Schwinger-boson Fock 
space as \<j,t) = a) a b\ |vac). Here and b\ create a state 
with spin and orbital directions a and r, respectively. 
Using these Schwinger bosons, we find that the state |^+) 
in Eq. (|4^|) can be written as 

N/2 

= n(4i-iAi-4-iAt) 



i=l 

x (a. 



2i,T U 2i+l,| u 2i,l u 2i+l,p 



|vac) , (52) 



and the interchange of a* and yields | V E'-). These 
staggered spin-orbital crystal states are first found by 
Kolezhuk and Mikeska^i The picture of these states is 
displayed in Figs. HJa) and|Ub). Obviously, the two- 
point correlation functions are nonvanishing only be- 
tween nearest-neighbor sites. For this spin-orbital 50(4) 
system, the string order parameters in Eq. (|47p can be 
written as 



km ((S?±T?)f[e i «W* r n(S*±T?)) = l (53) 



where (J, = x, y, or z. In the studies of two-leg spin 
ladders, these types of string order parameters were in- 
troduced to divide the topologically distinct gapped spin 
liquid states^ 

The fact that such two dimerized states are exact 
ground states of the projector Hamiltonian (fTTj) can be 
easily visualized when writing the projectors of the three 



50(4) channels in Eq. JlO]) as 

V k {i,j) = P S =o(i,j)Fk=o(i,j), 

= Ps=o(U)^t=i(m)+Ps=i(m)Pt=o(U), 
Vg(i,j) = P s=1 (i,j)P T ^(iJ), (54) 

where P s=0 (i,j) = ± - Sj • Sj and P s=1 (i,j) = f + 
Si ■ Sj are bond total spin projectors. Once the spin and 
orbital singlets are formed between nearest-neighbor sites 
in a staggered pattern, the 50(4)-symmetric projector 
Hamiltonian, 



H 



SO(4) 



Vg(i,i+1) 



= ^(S, • S i+1 + ^)(T, • T i+1 + 1), (55) 

i 

always annihilate such a spin-orbital crystal state. 

V. SO(n) BILINEAR-BIQUADRATIC MODEL 

As already mentioned, Hgo(n) is a bilinear-biquadratic 
Hamiltonian in terms of the SO(n) generators. More 
generally, we can also introduce a one-parameter family 
of the SO(n) symmetric bilinear-biquadratic model, 



Hi 



bb 



E 



cos e £ LfLllx + sin 6 ( £ if 

a<b \a<b j 



. (56) 

which is an extension of the familiar spin-1 bilinear- 
biquadratic model. The absence of the higher-order 
terms follows from the fact that such terms can be 
expressed via the lower-order terms by means of Eq. 
(fTU|) . To sketch the properties of this bilinear-biquadratic 
model, we need to identify several special integrable 
points. Let us introduce a slave boson representation, 



dl 



(57) 



which yields a constraint 53a=i dad a = 1. Using the slave 
bosons, the 50(n) generators can be written as L ab — 
i(d^ a dt, — d\d a ), and the 50(n)-singlet bond projector in 



Eq. pOjl is given by 



b- 



(58) 



Additionally, the SU (n)-invariant permutation operator 
is expressed as 



Q(iJ) = J2 d ^bd] b d ja - 



(59) 



ab 



Using the permutation operator and the singlet projec- 
tor, we can express the bilinear-biquadratic Hamiltonian 
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(156)) as 

H bb = ^{cos 9Q{i, i + 1) 

i 

+n[(n - 2) sin - cos 9]Pi(i, i + 1)}, (60) 

up to a constant. Several special points can be identified 
as follows: 

(1) 9 = tan -1 ^2 an d tan -1 — 7r; the Hamiltonian 
(|56p reduces to a sum of nearest-neighbor permutation 
operators, thus, it has an enhanced SU(n) symmetry. In 
this case, the transformation on each lattice site is in the 
SU(n) fundamental representation. For 9 = tan -1 — zq, 
this is the Uimin-Lai-Sutherland (ULS) model^ which 
can be solved by Bethe-ansatz method. It is known that 
there are gapless excitations above the ground states and 
the effective low-energy field theory is described by an 
SU{n)i Wess-Zumino-Witten models 

(2) 9 = the Hamiltonian (|5"6")l reduces to a sum 
of nearest-neighbor singlet projectors and also has an 
SU(n) symmetry. However, the transformations are in 
the SU(n) fundamental and its conjugate representations 
on the even and odd numbers of lattice sites, respec- 
tively. For 9 = — ^, a mapping to the n 2 -state quantum 
Potts model allows the model to be solved exactly ) 34 ' 35 i 36 
and the ground states are dimerized states with a finite- 
energy gap. 

(3) 9 = tan _1 (n - 4)/(n - 2) 2 ; the Hamiltonian (1561) 
was exactly solved by Reshetikhin21 via quantum inverse 
scattering method, which also exhibited gapless excita- 
tions. For n — 3, this point corresponds to the spin- 
1 Takhatajan-Babujian model^ which is the quantum 
critical point between Haldane gap phase and dimerized 
phase. For n = 4, the Reshetikhin point yields the 5(3(4) 
Heisenberg model, which is equivalent to two decoupled 
spin-1/2 Heisenberg antiferromagnetic spin chains. 

(4) 9 = tan -1 the ground states of the model Hamil- 
tonian ()56|) are just the matrix product states considered 
in Sees. Ill and IV. For an odd n, the ground state is a 
unique Haldane liquid state. For an even n, the ground 
states are twofold-degenerate dimerized states and are 
referred to non-Haldane liquid states. 

Therefore, these rigorous results suggest that an energy 
gap develops for the model (|56f in the finite parameter 
region 

tan- 1 < 9 < tan" 1 ^, (61) 

which always includes our matrix product ground-state 
point 9 = tan -1 —. The gap formation in this region is 
quite subtle. In the point of view of conformal field the- 
ory, the SU{n) ULS and the SO(n) Reshetikhin points 
are both conformal invariant and are characterized by 
two effective-field theories with different central charges. 
If so, there will be no renormalization flow from the 
SO(n) Reshetikhin point to the S'J7(n)-symmetric point 
according to Zamolodchikov's c theorem^ and an energy 




FIG. 5: The ground-state phase diagrams of the (a) 50(3), 
(b) 5*0(4), and (c) 5*0(5) symmetric bilinear-biquadratic spin 
chains. 



gap must be generated between these two conformal in- 
variant points. It was known that the conformal field the- 
ory for the ULS point is an SU(n)i Wess-Zumino-Witten 
model with central charge c = n — 1. The conformal field 
theory description for the SO(n) Reshetikhin point is an 
50(7j)i Wess-Zumino-Witten model with central charge 
c = n/2. In particular, the SO(3) Takhatajan-Babujian 
model is known to have a central charge c = 3/2 and the 
50(4) Heisenberg model (two decoupled spin-1/2 Heisen- 
berg antiferromagnetic spin chains) has a central charge 
c=2. 

Toward the odd n = 21 + 1 case, Itoi and Kato-- found 
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that a marginally relevant perturbation around the ULS 
point develops a Haldane gap for 9 < tan" 1 — z^, while 
the region 9 > tan -1 near the ULS point is mass- 
less. However, the even n case was extensively studied for 
n = 4, which corresponds to an 50(4) spin-orbital cou- 
pled systemi 41 i 42 i 43 i 44 i 45 These results reveal that there 
is a dimerized non-Haldane liquid phase with an en- 
ergy gap between the 5L/(4)-symmetric point and the 
50(4) Heisenberg point. In the non-Haldane liquid state, 
magnon excitations are incoherent and the low-energy 
excitations are a pair of solitons connecting two spon- 
taneously dimerized ground states. It is thus expected 
that these 50 (n) symmetric models also show such an 
interesting even-odd effect not only in the ground states 
but also in the low-energy excitations. For n = 21 + 1, 
the system is in a Haldane gap liquid phase with magnon 
excitations. For n — 21, the elementary excitations are 
solitons connecting the degenerate ground states. 

Following the exact results, the main phase diagrams 
of SO(n) bilinear-biquadratic model for n — 3, 4, 5 are 
displayed in Fig. ([5]). However, the SO(n) antiferro- 
magnetic Heisenberg model deserves more attention, cor- 
responding to the bilinear-biquadratic model (j56")l with 
a pure bilinear interaction for 9 = 0. When n = 3, 
it is just the quantum spin-1 antifcrromagnetic Heisen- 
berg model, which is in the Haldane gap region. When 
n = 4, the 50(4) Heisenberg model is equivalent to two 
decoupled spin- 1/2 antiferromagnetic chains, which have 
unique disordered ground states with power-law decay 
spin correlations. However, when n = 5, we find that 
the 50(5) antiferromagnetic Heisenberg model is not in- 
cluded in the Haldane gap region. Therefore, it is inter- 
esting to ask what are the ground states of the SO(n) 
antiferromagnetic models for n > 5. Based on a gener- 
alized Lieb-Schultz-Mattis theorem, Li 4 ^ studied SO(n) 
antiferromagnetic models for n = 4,5,6. He found that 
the 50(4) Heisenberg model is gapless, while 50(5) and 
50(6) Heisenberg models are suspected to have a gap. 
Together with our results, we predict that the SO(n) 
Heisenberg model for n > 5 belongs to the dimerized 
phase with a finite-energy gap. 



VI. CONCLUSION 

In conclusion, we have introduced a class of 
50 (n) symmetric spin chain Hamiltonians with nearest- 
neighbor interactions, whose exact ground states are two 
different SO(n) symmetric matrix product states de- 
pending on the parity of n. 

For an odd n = 21 + 1, a periodic chain has a unique 
ground state, which preserves an SO(2l + 1) rotational 
and translational symmetries. The 50(2/ + 1) symmetric 
spin chains with different I are directly related to quan- 
tum integer-spin chains belonging to the Haldane gap 
phase with a hidden antiferromagnetic order character- 
ized by nonlocal string order parameters. The hidden 
(Z2 x ^2)' symmetry responsible for the hidden order has 
been found by applying a unitary transformation to the 
model Hamiltonian. The Haldane gap and 4 l degenerate 
ground states in an open chain are natural consequences 
of this hidden symmetry breaking. 

For an even n — 21, a periodic chain has a twofold- 
degenerate dimerized ground state, which preserves 
50(2/) symmetry but breaks translational symmetry. 
These SO (21) matrix product states with different I are 
non-Haldane liquid states, which have soliton excitations 
connecting the two degenerate ground states. However, 
these SO(2l) matrix product states also contain a hidden 
antifcrromagnetic order characterized by nonlocal string 
order parameters. 

Finally, a generalized SO(n) symmetric bilinear- 
biquadratic model family has been discussed and the 
ground-state phase diagrams are sketched based on some 
known exact results. One of the important conclusions 
is that the ground state of the SO(n) symmetric Heisen- 
berg antifcrromagnetic spin model for n > 5 is predicted 
to be in a twofold-degenerate dimerized state. Further 
investigations on this are certainly required. 

We acknowledge the support of NSF of China and the 
National Program for Basic Research of MOST-China. 
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